Correlated dynamics of the motion of proton-hole wave-packets in a photoionized 

water cluster 



o 

(N 
O 

Q 

On 



Zheng Li (^ff), 1 Mohamed El-Amine Madjet, 1 Oriol Vendrell/'B and Robin Santra 1 ' 2 

1 Center for Free-Electron Laser Science, DESY, Notkestrasse 85, D-22607 Hamburg, Germany 
2 Department of Physics, University of Hamburg, D-20355 Hamburg, Germany 
(Dated: December 20, 2012) 

We explore the correlated dynamics of an electron-hole and a proton after ionization of a proto- 
nated water cluster by extreme ultra-violet (XUV) light. An ultrafast decay mechanism is found 
in which the proton-hole dynamics after the ionization are driven by electrostatic repulsion and 
involve a strong coupling between the nuclear and electronic degrees of freedom. We describe the 
system by a quantum-dynamical approach and show that non-adiabatic effects are a key element of 
the mechanism by which electron and proton repel each other and become localized at opposite sides 
of the cluster. Based on the generality of the decay mechanism, similar effects may be expected for 
other ionized systems featuring hydrogen bonds. 

PACS numbers: 82.20.Gk, 31.50.Gh, 33.80.Eh, 82.33.Fg 
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Correlated motions of electrons and protons in the 
so called proton coupled electron transfer (PCET) re- 
actions play a crucial role in charge transfer processes in 
a range of biological, chemical and material science sce- 
narios, and have been the subject of intensive study over 
many years [H-Q • These reactions can be thermally acti- 
vated processes occurring on relatively long time scales, 
in which the reorganization of the molecular environment 
plays a decisive role [HQ, as well as ultrafast PCET pro- 
cesses occurring upon photoexcitation Some exam- 
ples include the relaxation of electron vacancies in DNA 
and in crystals made of DNA bases, which are known 
to proceed via correlated proton and electron transfer 
reactions occurring in ps time-scales [5, 6]. Proton trans- 
fer (PT) processes triggered by a sudden change in elec- 
tronic structure are also known to protect DNA base- 
pairs from photo-damage caused by UV- visible light ab- 
sorption and mediate the decay back to the ground elec- 
tronic state 0- Non-adiabatic effects, in which the nu- 
clear motion strongly couples to the evolution of the elec- 
trons, are ubiquitous in the kind of processes described 
above. 

The response of matter to ionizing extreme ultravi- 
olet (XUV) radiation in the range 20 to 100 eV has 
been the subject of much interest in the last years due 
to opportunities offered by the generation of very short 
pulses in this energy range, opening the door to time- 
resolved studies of the induced ultrafast dynamics in 
the fs time regime The ionization from outer- 

valence shells can induce an ultrafast rearrangement of 
the electronic cloud in what are referred to as hole mi- 
gration processes [3, 17 1. The increased total charge 



and the modified bonding structure of a sample may re- 
sult in changes in its nuclear configuration reaching pos- 
sibly structures of strong non-adiabatic mixing between 
electronic states. Examples of ultrafast nuclear dynam- 
ics involving hydrogen atoms after extreme ultraviolet 
(XUV) photoionization are known and have been mea- 



sured by pump-probe techniques at the Free Electron 
Laser at Hamburg (FLASH) 8]. The fragmentation of 
protonated water clusters H(H 2 0)+ after XUV irradi- 
ation and removal of an outer-valence electron was re- 
cently measured at FLASH. The fragmentation products 
and their kinetic energies were recorded [Lsl . [l9j ]. Sys- 
tems in which a proton, and thereby an extra positive 
charge, is supported by a number of hydrogen-bonded 
molecules,are ubiquitous in soft condensed matter such 
as water [20( or in proton-exchan ge m embranes such as 
those used in fuel cell applications [2l| , where they medi- 
ate charge transport. A time-resolved picture and most 
mechanistic details of the nuclear-electronic dynamics af- 
ter their ionization, are however not yet known. 

Here, we investigate the decay and relaxation mech- 
anism of the electron-hole generated by absorption of 
XUV radiation in a prototypical water cluster, the so 
called Zundel cation, H(H20)J. We find an ultrafast 
coupled proton-hole relaxation process operating after 
photoionization from the outer- valence shell. The corre- 
lation function between the location of proton and hole 
shows that they respond to each other within few fs and 
quickly separate irreversibly. This is a consequence of 
the strong electrostatic repulsion between them, which 
results in potential energy surfaces (PES) for the nuclear 
motion with large gradients, but this is not the whole pic- 
ture. As discussed below, strong non-adiabatic effects in- 
volving the breakdown of the Born-Oppenheimer approx- 
imation play a crucial role in the relaxation mechanism. 
The fact that PES topologies are to a large extent related 
to the electrostatic repulsion between the hole and the 
proton indicates a general class of vibrational-electronic 
decay mechanisms operating in ionized hydrogen-bonded 
systems that should thereby be quite insensitive to par- 
ticular details of their molecular and electronic structure. 

The PES of the three lowest-energy states of the dica- 
tion are shown in Fig. QJi. They were calculated at the 
complete active space self-consistent field level of the- 
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FIG. 1. (color online) (a) Potential energy surfaces of the three lowest electronic states of H(H20) J + for the z (central proton 
position projected onto the oxygen-oxygen axis) and R (oxygen-oxygen distance) coordinates and keeping the rest of vibrational 
modes frozen. The dot marks the Frank-Condon nuclear geometry. The orbitals shown represent the highest occupied molecular 
orbitals at the Hartree-Fock mean field level, and they correspond to the electron-hole in the leading configuration for each of the 
three lowest energy electronic states of the dication in a multiconfigurational description, (b) Scheme of the hydrogen-bonded 
system indicating the location of the electron pairs that can be ionized and the corresponding potential energy surfaces as a 
function of the central proton position. Ionization leading to the electronic states under consideration can take place from lone 
electron pairs on each monomer (in red and blue) or from the electrons involved in the hydrogen bonding (in green). 



ory [22J with the Molcas suite of programs [23|. The 
two nuclear coordinates of the 2D cut in Fig. QJi cor- 
respond to the oxygen-oxygen distance (R) and to the 
projection of the central proton position onto the oxygen- 
oxygen axis (z), whereas the rest of coordinates are held 
fixed. A schematic representation of the vibrational co- 
ordinates may be found in Fig. [2J The point in coordi- 
nate space marked with a dot signals the position of the 
minimum energy configuration in the ground electronic 
state of H(H20)J. Therefore, the nuclear wavepacket 
is centered at this position immediately after photoion- 
ization. At such geometry of the cluster the two low- 
est electronic states of the dication are energetically de- 
generate [24]] . The orbitals shown in Fig. QJi indicate 
the electron-hole in a 1-particle picture. The electronic 
configurations with the missing electron in such orbitals 
contribute to more than 90% in the respective multicon- 
figurational electronic wave functions, indicating that a 
1-particle picture description of the states under consid- 
eration is adequate in this case. By inspecting them we 
realize that the two lowest-energy states can be described 
as removing an electron from a lone electron pair in cither 
water molecule. 

The topology of the three low lying PES of the dica- 
tion is illustrated schematically in Fig. [T]d as a function 
of the central proton position. The shape of the PES is 
mostly the result of the electrostatic repulsion between 
the positive charge density associated with the central 
proton and that of the electron-hole. This is made clear 
by noting that the gradient of the PES at the central 
point of conical intersection is about 6.5 eV-A -1 . This is 
equivalent to two point charges of equal sign at a distance 
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FIG. 2. (color online). Schematic representation of the 7 
internal nuclear coordinates used in the wavepacket propa- 
gations showing in each case the relevant atoms. The pro- 
ton transfer coordinate used in this work was transformed as 
z' — z/(R — S) ,25j], where 8 is a constant. Along this work 
appearances of (z) refer to the (z') defined here. 



of 1.4 A; the distance of the central proton to one of the 
oxygen atoms is about 1.2 A at the equilibrium geome- 
try of the cluster. When the central proton is displaced 
towards e.g. the water molecule on the left, the elec- 
tronic state with the hole on the water molecule on the 
right side stabilizes, while the state with the hole on the 
left side destabilizes, resulting in the conical intersection 
between the electronic states as a function of the proton 
position. The third electronic state of the dication can be 
described by an electron-hole delocalized mostly around 
the central proton. Therefore, the PES is quite flat along 
the PT coordinate in the vicinity of the Franck-Condon 
point. When the proton is displaced towards either water 
molecule, however, the electron-hole eventually becomes 
more stable on a localized water site, and at that geome- 
try the third dicationic state crosses with the second one, 
resulting also in a seam of conical intersection. 

The ultrafast dynamics of the protonated water dimer 
unfolding after outer-valence ionization were computed 
quantum dynamically using the Heidelberg implemen- 
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tation [26| of the multiconfiguration time-dependent 
Hartree method [27L[28j|. Out of the total 3N a - 6 = 15 
vibrational coordinates of the system, where N a equals 
the number of atoms in the system, seven were kept ac- 
tive for the dynamics calculations and the other eight 
were kept frozen at their ground state equilibrium val- 
ues. The active vibrational modes are shown in Fig. [2] 
and correspond to the three Cartesian coordinates lo- 
cating the central proton in the frame defined by the 
water molecules, the oxygen-oxygen distance, the wa- 
ter molecule wagging modes, and the relative angle be- 
tween both water molecules around the axis connecting 
both oxygen atoms. The remaining frozen coordinates 
correspond to the internal degrees of freedom of each 
monomer, which are considered as rigid bodies, and the 
two rocking angles. The separation in active and inac- 
tive modes is made based on the vibrational frequencies 
of the involved motions. Higher frequency modes are kept 
frozen, since they are expected to be less coupled to the 
geometric rearrangements after ionization. A detailed 
description of the coordinates used and of the kinetic en- 
ergy operator in this set of internal coordinates may be 
found in Ref. @. 

The Hamiltonian of the system in diabatic representa- 
tion reads 12911 



H = f n ®l d + ^\i)Wi,j(Q,){j\ 



(1) 
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where Wij operates on the space of the vibrational de- 
grees of freedom, Id = J2i i s the unit operator in 
the discrete space of the diabatic electronic states, and 
|z) refers to the ith diabatic electronic state. The diaba- 
tization procedure used is based on regularized diabatic 
states [30j and details will be presented elsewhere. Af- 
ter ionization, assuming a sudden electron ejection that 
leaves the system in the ith electronic state, the corre- 
sponding wavepacket is ^(i = 0) = $ (Q)®|i). The nu- 
clear part $o(Q) corresponds to the ground vibrational 
state on the ground electronic state PES of H(H20)J. 
Such wavepackets are propagated for every initial elec- 
tronic state in coupled PES using the Hamiltonian in 
Eq. |T]). We assume the ionization cross-section to be 
constant over the considered energy range for the initial 
electronic states. Therefore, the observables calculated 
from each propagation are averaged with equal weight. 

The ultrafast correlated dynamics of hole and proton 
are made explicit by inspecting the correlation function 
between their positions, shown in Fig. [3] The correla- 
tion function is defined as C p h = S p h / \J 'S pp Shh, with 
S X y = (xy) — {x){y). C p h is bound to take values between 
— 1, complete anti-correlation, and +1, complete correla- 
tion. For the excess proton position operator we simply 
take z, the projection of the central proton position onto 
the oxygen-oxygen axis of the cluster. A convenient def- 
inition of the hole position operator that can be readily 
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FIG. 3. Correlation function between proton and electron- 
hole positions obtained from the propagated wavepacket. 



applied to the propagated wavepacket is 
h = ±R®(\r){r\-\l){l\), 



(2) 



where \l) and |r) refer to the two diabatic electronic states 
with the hole localized on the left- and on the right-hand 
side water monomers, respectively, and R is the posi- 
tion operator corresponding to the oxygen-oxygen dis- 
tance. According to this definition, the average location 
of the electron-hole corresponds to the position of the left 
(right) oxygen atom when the system is in the electronic 
state \l) (|r)). The correlation function in Fig. [3] is ob- 
tained by summing over the different initially populated 
electronic states. At t = the proton and hole positions 
are still uncorrelated. The degree of anti-correlation in- 
creases very quickly as the proton is displaced towards 
one of the water monomers and the hole moves in the 
opposite direction, reaching almost —1 in about 5 fs. A 
small recurrence in the correlation function is seen be- 
tween 7 and 11 fs. It originates from the bounce-back 
motion of the proton against the water monomers, which 
have practically not yet started moving simply due to 
their higher inertia. After 10 to 20 fs proton and hole 
have irreversibly localized at opposite sides of the sys- 
tem, although the average oxygen-oxygen distance has 
only increased by 0.2 A and larger recurrences could be 
expected. This irreversibility is a direct consequence of 
the multidimensional nature of the triggered vibrational- 
electronic motion, which prevents the wavepacket, just 
a few fs after photoexcitation, from returning to nuclear 
and electronic configurations resembling the initially un- 
correlated state. 

The ultrafast nuclear dynamics and fragmentation af- 
ter ionization are shown in Figs. H^,, HJd and Ht, which 
illustrate the dynamics of the cluster in the z, R coordi- 
nates after integration over all other nuclear coordinates 
and electronic states. As the proton and hole quickly 
separate, the central proton approaches the vicinity of 
an oxygen atom in less than 5 fs. At about 10 fs there is 
a revival in which the proton bounces back to the center 
of the cluster because, as already mentioned, the oxygen 
atoms did not yet have time to move. After 15 fs the oxy- 
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FIG. 4. (color online). Reduced probability densities in the z (central proton position projected onto the oxygen-oxygen axis) 
and R (oxygen-oxygen distance) coordinates obtained from the propagated wavepackets after tracing over all other vibrational 
modes and electronic states. Snapshots (a), (b), and (c) correspond to propagations with the full Hamiltonian in Eq. JTJ. 
Snapshots (d), (e), and (f) correspond to propagation without non-adiabatic coupling in the Hamiltonian, Eq. Q. 



gen atoms have already separated by about 0.2 A, and 
the probability that the proton is found at the central 
position between the water molecules has almost van- 
ished. At this point, the correlated dynamics between 
the proton and the electron-hole are over and they are 
found localized at opposite sites of the cluster, leading to 
fragmentation of the system with 100% probability. 

In order to illustrate the dramatic role of the non- 
adiabatic coupling (NAC) between nuclear motion and 
electron-hole dynamics, we briefly discuss the quantum 
dynamics of the dication upon ionization in the adia- 
batic representation, but neglecting the NAC altogether. 
Therefore the wavepackets evolve in uncoupled Born- 
Oppenheimer PES of the ionized system and the uncou- 
pled Hamiltonian operator reads 

H = f n ®i a + J2\a)V a (Q)H (3) 

a 

where l a = ^ Q |a)(o;| and |a) are now uncoupled adi- 
abatic electronic states. The reduced probability den- 
sity p(z, R) in which all other nuclear coordinates and 
electronic states are integrated over is shown for this 
case in Figures [4ji, Hk, and Hf. The part of the ini- 
tial wavepacket on the ground electronic state PES of 
the dication, i.e. initially below the lowest energy coni- 
cal intersection, escapes from the Franck-Condon region 
in about 5 to 10 fs, quickly leading to the dissociation 
channel H30 + + H2O" 1 ". The reason is that in the low- 
est adiabatic electronic state of the dication, the central 
proton and electron-hole are always at opposite sides of 
the cluster (cf. Fig. [IJ. This is not the case when the 
nuclear dynamics start in the two lowest excited states 
of the dication. Although the cluster is doubly charged, 
the missing coupling between nuclei and electrons pre- 
vents it from fragmenting, the electron-hole and proton 
cannot rearrange and the system stays bound. This is 
clearly seen by comparing Figs. 2k and Hfl. It is worth 



mentioning at this point that the fragmentation exper- 
iments in Ref. [18] do not report unfragmented clusters 
after ionization. 

In conclusion, we have described an ultrafast decay 
and relaxation mechanism by which an outer- valence hole 
generated by XUV light in a strongly hydrogen-bonded 
system triggers proton-hole correlated dynamics unfold- 
ing in a time-scale of just a few fs. The motion of the 
proton and the hole is anticorrelated and strongly cou- 
pled and their dynamics are driven by their electrostatic 
repulsion, which results in PESs with large relative gra- 
dients at regions of intersection and strong non-adiabatic 
effects. The rearrangement of the proton and hole leads 
to electronic relaxation of the system, indicating that the 
dynamics of the lightest nuclei is important when consid- 
ering the motion of electron-holes in molecular systems 
and clusters. Artificially removing the non-adiabatic cou- 
pling in our quantum dynamical calculations completely 
quenches the described mechanism because proton and 
hole cannot rearrange. The electrostatic nature of the 
reported mechanism suggests that this class of decay 
pathway may be a common feature of ionized hydrogen- 
bonded systems. 
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